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Abstract 

In the framework of a continuum theory, it is shown that the direct flexoelectric response of a 
finite sample essentially depends on the surface polarization energy, even in the thermodynamic 
limit where the body size tends to infinity. It is found that a modification of the surface energy can 
lead to a change of the polarization response by a factor of two. The origin of the effect is an electric 
field produced by surface dipoles induced by the strain gradient. The unexpected sensitivity of the 
polarization response to the surface energy in the thermodynamic limit is conditioned by the fact 
that the moments of the surface dipoles may scale as the body size. 
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The impact of the sample surface on the macroscopic properties of ferroics is an issue 
that has been attracting the attention of workers during the past 50 years. Starting from the 
seminal papers for magnetic [l| and non-magnetic ferroics, which were based on Landau 
theory, this issue latter became a matter of extensive hrst principles studies (see e.g. jj]). 
All obtained results suggest a pronounced size effect, i.e. an increasing of sample surface 
impact on its macroscopic properties while decreasing the sample size. Quantitatively, strong 
effects were predicted for samples having small dimension not exceeding tens of nanometers. 

A good example of such phenomena is the out-of-plane dielectric response of ferroelectric 
thin hlms. Here, depending on the film thickness, the effective dielectric constant acquires 
a relative correction of the order oi C/h, where h is the film thickness and C is the microscopic 
scale, which typically does not exceed a few tens of nanometers (see review paper jo] and 
references therein). 

An essential feature of such size effect is that its strength is, as expected, sensitive to 
the properties of the surface of the sample, specifically to the value of the polarization 
dependent contribution to the surface energy. For example, if such a contribution vanishes 
(the so-called free boundary conditions) then the aforementioned correction vanishes as well. 
In the thermodynamic limit h —)■ oo this correction vanishes also. 

The trends in the size effect outline above are of interest not only from the point of view 
of fundamental science, they also provide guidelines for down-scaling of microelectronic and 
micro-mechanical devices. In this context, the verihcation of these trends for electrome¬ 
chanical effects in solids seems to be a problem deserving a theoretical treatment. In this 
letter we present a simple modelling of the impact of the surface polarization energy on the 
average polarization response of a plate to its bending, which is associated with the bulk 
flexoelectric effect. We show that, in contrast to the dielectric response, a modihcation of 
the surface polarization energy can readily lead to a 100 % modification of the flexoelectric 
response of a hnite sample, which holds in the thermodynamic limit h ^ oo . Particularly, 
we demonstrate that, at h —)■ oo, a modihcation of the surface polarization energy can lead 
to a change of the polarization response by a factor of 2. 

The hexoelectric effect, as a polarization response to a strain gradient, has been originally 


introduced as a purely bulk effect being viewed as a simple analogue of the piezoelectric 
effect j?, ^ 
the problem 0 


or a historical overview see e.g.j9|). However, further theoretieal studies of 
revealed the presence of unexpected contributions having no analogues 
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in piezoelectricity. It was, for instance, realized that the flexoelectric response of a hnite 
sample should be essentially influenced by effects controlled by the sample surface. This 
result, originally based on general arguments, recently received a direct conhrmation based 
on hrst principles calculations [n|. 

However, in the phenomenological framework one can single out the so-called static bulk 
flexoelectric effect, which can be viewed as an analogue of the piezoelectric effect and others, 
additional contributions {q]. It is this effect that is customarily used for the description of 
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the flexoelectricity-related phenomena in continuum theory of hnite samples (see e.g. 

4). It is usually addressed in terms of the following density of a thermodynamic potential 

Im\ 


i29ijklPi,j Pk,l T T 2 , (1) 


where Ei, Pi, and Uij are the Cartesian components of electric held, polarization and strain 
tensor respectively, the suffixes separated by a comma mean the corresponding spatial deriva¬ 
tives. The hexoelectric ehect is controlled here by a Lifshitz-type invariant containing the 
so-called hexocoupling tensor fijki- Hereafter the Einstein summation convention is adopted. 

For the description of the direct hexoelectric response of a hnite sample, i.e. polarization 
response to a strain gradient, similar to the case of the dielectric response (c.f. j^), one 
should minimize with respect to the polarization the total thermodynamic potential of the 
sample, which incorporates the surface polarization energy introduced by Kretschmer and 
Binder [^: 

F = / TvdV + / \alPiPjdS. (2) 


^ -flflPidS. 

Here, in the last rhs term, the integration is done over the sample surface and the two-suffix 
quantity afp controlling the surface polarization energy, is dependent on the orientation of 
the surface. 

Minimization F with respect to the polarization yields the equation of state 


^IkPl 9ijklPi,j,l fijkl^ij,l Ek 0 

with the boundary conditions QQ 

^IkEl ffi 9ijklPi,j^l T 'zfijkl'^ij'kf'l 0 , 


( 3 ) 


( 4 ) 


where rii are the components of the unit vector normal to the surface. 
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The framework presented above enables the description of both dielectric and flexoelectric 
responses of a hnite sample. For the geometry of the parallel plate capacitor the former was 
addressed in many papers (see e.g. [j, la, ll^). A remarkable qualitative conclusion of this 
theoretical analysis is the difference between the dielectric response in the limiting cases 
where = 0 (the so-called free boundary condition) and ^ oo (blocking boundary 
condition). In the ’’free” case, where the polarization surface energy vanishes, the effective 
value of the dielectric constant was found to be independent of the plate thickness h and 
equal to its bulk value. At the same time, in the ’’blocking” case the effective value of the 
dielectric constant was found to be smaller than its bulk value, with the difference scaling 
as l/h. An important feature of these results is that in the ’’thermodynamic limit” h ^ oo 
the dielectric response becomes insensitive to the polarization surface energy. 

In the present letter we test the sensitivity of the flexoelectric response of a plate, con¬ 
ditioned by the static bulk contribution, to the polarization surface energy. To do this, we 
calculate it for the free and blocking cases to hnd a situation drastically different from that 
of the dielectric response. Specihcally, we hnd that (i) in the free case the response can be 
weaker that in the blocking one, (ii) the difference between the cases does not disappear in 
the thermodynamic limit h ^ oo and (iii) in the limit h ^ oo the blocking case yields the 
hexoelectric response twice that of the free case. 



FIG. 1. A central cross-sections of plate: dashed line - before bending and solid line - after. 


Consider a thin (001) plate of a cento-symmetric cubic material of thickness h normal 
the X 3 Cartesian coordinate, the origin being placed in its center (Fig. [T]). The plate is 


symmetrically bent, so that the curvature of the p 
normal to it. Following a thin plate approximation 

C1122 


ate G is the same in all cross-sections 


Mil — U22 — X^G; M33 — —2X3 


Cllll 


17l | we present the stains in the form 
G] U 12 = U 23 = Mi3 = 0, (5) 
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where C1122 and cim are the components of the tensor of elastic moduli. In the same 
approximation, in view of the geometry of the system and the symmetry of the material, all 
variables depend only on the Cartesian coordinate and only P3 and P3 components of 
the polarization and electric held, denoted hereafter as P and E, respectively, are non-zero. 

Using Eqs. (E]) and (E]) and denoting the derivative with respect to X3 as the prime, the 
polarization equation of state reads 

F = Fg + xgP" + XE , (6) 


where g = 5(3333, y = 1/033, and 


Ciiii/1122 — Cii22/llll„^ 

Fg = X- 2 G. 

Cllll 


(7) 


It is clear from Eq. ([6]) that Pq signihes the polarization, induced by the plate bending, if 
the polarization were homogeneous and the electric held in it were zero. 

For the general situation one should take into account the Poisson equation for the elec¬ 


trical displacement D, which we presen 
depolarizing ehects in ferroelectrics jo. 


in the form allowing an adequate description of the 
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D — P e^E, 


( 8 ) 


where £b is the background permittivity while F is redehned as the ” hexoelectric” contribu¬ 
tion to the polarization. 

Since all variables are functions of the coordinate X3 only, the Poisson equation implies 
D = const. So that 

(F) -|- £h{E) = P + EhE, ( 9 ) 


where 

(E) 

Combining Eqs. (jb]) and 


h/2 

^ j EdX^ and (F) 

-hl2 

we hnd 


hl2 

i / PdX,. 

-hl2 


( 10 ) 


(1 + - XgP" = xhh{P) + Pg + x{E)- (11) 


Equation flTTll appended with the boundary conditions specify the problem. Using Eq. ( 01 ), 
the latter read 


F 


X3=±/i/2 


0 


( 12 ) 
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for the blocking case {aff. —)• cx)) and 


P'\x3=h/2 — —P'\x->=-h/9. — 


h 


X3=-hl2 


^9X 


G 


( 13 ) 


for the free case (a^ = 0). 

The solutions to linear differential equation flTT]) . satisfying boundary conditions (IT^ or 
uni), can be routinely obtained. However, being interested in (P) only and using the results 
available in the literature, one can evaluate the polarization response sought without explicit 
solving of the differential equation, as presented below. 

Consider the blocking case. The problem of the direct flexoelectric response of the shorted 
plate with the blocking boundary conditions, given by Eqs. (ITT]) and flT^ . can be readily 
solved using results on the dielectric response of the plate with such boundary conditions. In¬ 
deed, as is clear from Eq. fITT]) the flexoelectric response to the strain gradient, characterised 
by a value of Pq, is equivalent to the dielectric response to the average held (P) = Pg/x, 
calculated neglecting ferroelectricity. Thus, the average polarization produced by the strain 
gradient can be cast in the form 

{P) = Pg—, ( 14 ) 

X 


where Xes is the effective dielectric susceptibility of the sys 
neglecting hexoelectricity. Results of the earlier papers {2, 6 


em. 


16 | 


i.e. (P) / (E) is calculated 
suggest 


1 - 2^/h 


(15) 


where 




9 X 


( 16 ) 


1 + x/^h 

Equation (IT^ is valid under the condition h » ^. Using typical values of g in ferro- 
electrics, one hnds that, depending on the value of x/^b, ^ varies in the interval between the 
lattice constant and the correlation length of the material, which is about . Thus, this 
condition should hold for any realistic physical situation. Equations (ITTl) and (lT 5 l) implies 

1 - 2^/h 


(P) = Pc 


1 + (2f//i)(x'/£b)' 


(17) 


Note that (P) = Pq in the ’’thermodynamic limit” h —)■ cx3. 

In contrast to the blocking case, for the free case the problem cannot be reduced to that 
for the dielectric response of a plate having vanishing surface polarization energy since the 


6 










boundary conditions for the latter, P'\xi=±h/2 = 0 are different from flTSD . However, by 
averaging Eq. flTT]) over the plate thickness with {E) = 0 , we hnd a relationship 

(-P) = -Pg + ^ {P'\xi=h/2 - P'\xz=-h/2) , ( 18 ) 

which being combined with Eq. flT^ yields the flexoelectric response of the shorted plate as 


(P) = Pg/2. 


( 19 ) 


flifil- 


It is instructive to recall that in the free case Xes = X 

The above results reveal a drastic difference between the impact of the surface polariza¬ 
tion energy on the dielectric and flexoelectric responses. The most striking is the fact that 
in the thermodynamic limit h ^ oo the flexoelectric response is two times different between 
the blocking and free cases while there is no such difference for the dielectric response. The 
behaviour of the dielectric response is readily expected: the impact of a surface perturba¬ 
tion (polarization surface energy) vanishes in the thermodynamic limit. Meanwhile, even 
in the thermodynamic limit the contribution of the bulk flexoelectric effect to the direct 
flexoelectric response was found to be sensitive to the surface perturbation. Actually, a 
similar phenomenon was identihed in the case of the contribution of the so-called surface 
piezoelectricity to the total flexoelectric response of the sample 9 |, ll 9 (]. Like in the case of 
surface piezoelectricity, under short-circuited condition, the strain gradient induces surface 
dipoles which in turn create an electric held in the body of the plate. The formation of such 
dipoles is seen in Fig. [ 2 l The surface dipoles form in both blocking and free cases. However, 
in the former, the surface dipole moment is independent of the plate thickness h while in 
the latter it scales as h. In the free case, this gives rise to a non-vanishing impact of the 
surface perturbation in the thermodynamic limit. 

It is worth to mentioning that, in view of basic thermodynamic arguments, one expects 
that the two-times difference between the direct hexoelectric effects in the blocking and 
free cases must imply the same difference between the converse flexoelectric effects (bending 
of the plate in a homogeneous electric held) in these cases. This expectation is in full 
agreement with the calculations of the held-induced plate bending (converse hexoelectric 
ehect) controlled by the static bulk hexoelectric ehect performed for the blocking 19 | and 


free 

case. 


20 l | cases. Indeed, the converse ehect was found to be twice stronger in the blocking 
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FIG. 2. Schematics for the distribution of polarization across plate thickness, caused by its bending: 
(a) blocking boundary conditions, (b) free boundary conditions, h and ^ and the plate thickness 
and the spatial scale given by Eq. m, respectively. 


One should also comment on the applicability of the continuum approach to the prob¬ 
lem addressed. The typical spatial scale for the polarization variation is ^ as clear from 
Eqs. (HH), (USD, and Fig. [ 2 j The applicability of the continuum approach requires that ^ 
should appreciably exceed the lattice constant of the material. This condition is readily met 
for the so-called weak ferroelectrics [l^, 21], where is of the order of y and ^ is about the 
correlation length that, in turn, by dehnition appreciably exceeds the lattice constant. In 
the case of ’’regular” proper ferroelectrics, where £b << X) ^ is expected be about the lattice 
constant of the material. This pushes the continuum approach to the edge of its applicability 
range. However, one may expect that the qualitative conclusions of the presented analysis 
will hold in this case also. 

In conclusion, we have shown that the direct flexoelectric response of the plate associated 
with the static bulk flexoelectric effect is very sensitive to the surface polarization energy. 
In contrast to the polarization response, a modihcation of the surface polarization energy 
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can result a variation of the flexoelectric response by a factor of 2. Remarkably, the impact 
of such surface perturbation holds in the thermodynamic limit. The reported results sug¬ 
gest that any modelling of electromechanical hnite systems involving flexoelectricity should 
take into account the surface polarization energy and the modihcation of the polarization 
boundary conditions associated with the flexoelectric coupling. Nowadays such modellings 
customarily neglect these issues. 

This project was supported by the Swiss National Science Foundation. The Research was 
also supported by grant of the government of the Russian Federation 2012-220-03-434. 


[1] M. I. Kaganov and A. N. Omelyanc, Sov. Phys. JETP 34, 895 (1972). 

[2] R. Kretschmer and K. Binder, Phys. Rev. B, 20, 1065 (1979). 

[3] J. Junquera and P. Ghosez, Nature (London) 422, 506 (2003). 

[4] M. Stengel, D. Vanderbilt, and N. A. Spaldin, Nature Mater. 8, 392 (2009). 

[5] It is calculated from the derivative of the charge of the plates of a parallel plate capacitor 
containing the film with respect to the applied voltage. 

[6] A.K. Tagantsev, G.Gerra, J. Appl. Phys. 100, 051607 (2006). 

[7] S. M. Kogan, Sov. Phys. Solid State 5, 2069 (1964). 

[8] E. V. Bursian and N. N. Trunov, Sov. Phys. Solid State 16, 760 (1974). 

[9] P. V. Yudin and A. K. Tagantsev, Nanotechnology 24, 432001 (2013). 

[10] A.K. Tagantsev, Physical Review B, 34, 5883 (1986). 

[11] M. Stengel, Physical Review B, 90, 201112, (2014). 

[12] E.A. Eliseev, A.N. Morozovska, M.D. Glinchuk, R. Blinc, Phys. Rev. B 79, 165433 (2009). 

[13] M. S. Majdoub, P. Sharma, T. Gagin, Phys. Rev. B 77, 125424 (2008). 

[14] Z. Zhang, Z. Yan, and L. Jiang, J. Appl. Phys. 116, 014307 (2014) 

[15] These boundary conditions were originally derived in [l^ . where these are given with a mis¬ 
print in the order of suffices of the tensor gijki, corrected in the present paper. 

[16] O.G. Vendik, S.P. Zubko, J. Appl. Phys. 88, 5343 (2000). 

[17] L. D. Landau and E. M. Lifshitz, Theory of Elasticity (Pergamon Press, Oxford, 1975). 

[18] A. K. Tagantsev, Ferroelectrics, 375, 19 (2008). 

[19] A.K. Tagantsev, A.S. Yurkov, J. Appl. Phys. 112, 044103 (2012). 


9 



[20] A.S. Yurkov, Phys. Solid State 57, 460 (2015), (in press). 

[21] A. K. Tagantsev, I. G. Sinii, and S. D. Prokhorova, Izvestiya Akademii Nauk SSSR Seriya 
Fizicheskaya 51, 2082 (1987). 


10 



